Single polymer dynamics in elongational flow and the confluent Heun equation 
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We investigate the non-equilibrium dynamics of an isolated polymer in a stationary elongational 
flow. We compute the relaxation time to the steady-state configuration as a function of the Weis- 
senberg number. A strong increase of the relaxation time is found around the coil-stretch transition, 
which is attributed to the large number of polymer configurations. The relaxation dynamics of the 
polymer is solved analytically in terms of a central two-point connection problem for the singly 
confluent Ifeun equation. 
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I. INTRODUCTION 



Dilute polymer solutions exhibit physical behaviours that distinguish them from ordinary Newtonian fluids. Even 
a small polymer concentration can considerably change the large-scale behaviour of the flow by enhancing viscosity 
' and reducing the turbulent drag. A comprehensive understanding of hydrodynamical properties of polymer solutions 
is still lacking despite the large number of industrial applications (e.g. Sreenivasan and White 2000). A starting 
point for the theoretical description of dilute polymer solutions is the dynamics of isolated polymer molecules. The 
knowledge of how a single molecule is deformed by the velocity gradient allows the development of constitutive models 
I I - that in turn can be used to calculate the large scale flow. The stationary dynamics of a single, isolated polymer has 
received much attention to date (see Larson 2005, Shaqfeh 2005, for a review). On the contrary, less is known 
, about non-equilibrium dynamics of isolated polymers in flow. Since Rouse (1953) and Zimm's (1956) seminal works, 
■ experimental, numerical, and theoretical studies focused on the internal relaxation dynamics of a polymer floating 
[ in a solvent under the influence of Brownian fluctuations. The two model situations considered were: (a) a polymer 
• suspended in solution and pulled at the ends (Quake et al. 1997, Hatfield and Quake 1999); (b) a tethered polymer 
^\ ' submitted to a uniform flow and freely relaxing after cessation of the flow (Perkins et al. 1994, Brochard-Wyart 1995, 
, Manneville et al. 1996, Rzehak and Zimmermann 2002). The above studies did not take into account the interaction 
' between the polymer and an external flow. Hernandez Cifre and Garcia de la Torre (1999) note that Rouse and 
Zimm's theories may not provide the adequate time scale for coil-stretch processes in flow on the basis of Brownian 
(~| Dynamics simulations. Here, we investigate how an elongational flow influences polymer relaxation dynamics. We 
* ^ determine the time scale associated with polymer deformation in the flow and show a significant deviation from Rouse 
(Ih ' and Zimm's predictions in the vicinity of the coil-stretch transition. Our analysis is based on an eigenvalue problem 
^ I for the confiuent Heun equation, and constitutes a new physical application of this latter equation. 

Polymer dynamics is extremely rich already in simple deterministic flows, such as elongational flows (Perkins et 
al. 1997, Smith and Chu 1998), shear flows (Smith et al. 1999, Celani et al. 2005b, Puhafito and Turitsyn 2005, 
Gerashchenko and Steinberg 2006) or combinations of the two (Hur et al. 2002, Babcock et al. 2003). The velocity 
gradient of a non- uniform flow stretches the polymer, while entropic forces attempt to restore the polymer into the 
coiled equilibrium shape. In other words, the configuration of the polymer results from the counterbalance between 
the entropic forces and the hydrodynamic drag. Here we consider an elongational fiow, which is defined by a constant 
velocity gradient A. This flow is particularly effective in stretching polymers far from their equilibrium configuration. A 
transition to the stretched state of the polymer occurs as the velocity gradient exceeds the critical value Ac = 1/(2t), 
where r is the time associated with the slowest relaxation mode of the polymer in thermal equilibrium with the 
surrounding medium. For A < Ac polymers stay in the coiled equilibrium configuration; for A > Ac they become fully 
extended. This phenomenon is known as the coil-stretch transition (de Gennes 1974). Rouse (1953) computed r as a 
function of the number of monomers in the polymer assuming the polymer could be described as a beads-and-springs 
chain with Hookean interactions. Zimm (1956) refined Rouse's prediction by taking into account hydrodynamic 
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interactions between segments of the chain. 

The study of polymer dynamics has benefited from a class of mesoscopic models that are based on a coarse- 
grained treatment of the polymer molecule. One of the simplest is the elastic dumbbell model, which only takes into 
account the slowest oscillation mode of the molecule (Bird et al. 1987). Notwithstanding this crude simplification, 
the dumbbell model captures the main aspects of polymer dynamics in elongational fiows, such as the coil-stretch 
transition (Perkins et al. 1997) and finite-time conformation hysteresis (Schroeder et al. 2003, Schroeder et al. 2004). 
We examine how the probability density function (PDF) of the extension of a dumbbell-like molecule approaches its 
stationary form. We compute the typical time it takes for the initial PDF to attain its steady-state form and show 
that, in the proximity of the coil-stretch transition, this time is exceedingly long compared to r. A similar behaviour 
is encountered in white-in-time isotropic random flows (Celani et al. 2005a, Martins Afonso and Vincenzi 2005); 
however, for an elongational flow the amplification of the transient relaxation time is much stronger. 

The problem is solved within the framework of the Fokker-Planck equation for the PDF of the extension of the 
polymer. The computation of the transient relaxation time is recast as a central two-point connection problem (CTCP) 
for a singly confiuent Heun equation. The Heun equation is the general Fuchsian differential equation with four 
regular singularities. Its singly confiuent form, also known as the generalised spheroidal wave equation, results from 
the merging of two regular singularities into one irregular singularity of Poincare rank 1 (see Slavyanov 1995, for a 
review). A CTCP is an eigenvalue problem for an ordinary differential equation where: (a) at both endpoints of 
the interval of definition is located a singularity of the equation; (b) the solutions arc required to have a specified 
asymptotic behaviour while approaching the two singularities from inside. Contrary to the hypergeometric equation, 
an explicit formula for the CTCP is not known for the Heun family of equations. The eigenvalues can be determined 
only as solutions of transcendental equations involving continued fractions. Well-known applications of the confiuent 
Heun equation are the electronic spectra of the hydrogen-molecule ion in the Born-Oppcnheimer approximation (Jaffe 
1934) and the Teukolsky equation describing small perturbations to the Kerr geometry in black hole theory (Leaver 
1985). For a comprehensive introduction to physical applications of the Heun equations we refer the reader to the 
book by Slavyanov and Lay (2000). 

The rest of this paper is organised as follows. The dumbbell model is introduced in section [TTl The computation 
of the transient relaxation time and the related CTCP for the confluent Heun equation are presented in section IHII 
Section [TVl is devoted to conclusions. 



When one is interested in the statistics of polymer extension the dynamics of a polymer molecule can be described, 
to a first approximation, in terms of its slowest oscillation mode. In this case, the molecule can be modelled as an 
elastic dumbbell, i.e. as two beads connected by a spring. The beads represent the ends of the molecule and the spring 
models entropic forces. The separation between the beads measures the extension of the polymer. When introduced 
into a non-uniform fiow the molecule experiences collisions with fluid particles and becomes stretched under the 
action of the velocity gradient. In the simplest case, the drag force is assumed to be proportional to the velocity of 
the polymer relative to the fluid and thermal agitation is modelled by Brownian motion. In most applications the 
extension of the polymer remains smaller than the dissipative scale of the carrier fiow, and therefore the velocity 
field can be assumed to be linear^. The separation vector of the ends of the polymer, Q, is then a stochastic process 
evolving according to the three-dimensional stochastic ordinary differential equation^ (e.g. Bird et al. 1987) 



where Q = \Q\, Qq is the equilibrium mean-square separation of the ends of the molecule, and B{s) denotes the 
three-dimensional standard Brownian motion. The function / determines the entropic force. We consider the finitely 
extensible nonlinear clastic (FENE) dumbbell model with 



II. ELASTIC DUMBBELL MODEL 




(1) 



f{Q) 



1 



< g < L, 



1 - Q^/V 



2 ' 



^ Davoudi and Schumacher (2006) recently investigated the situation where the extension of the polymer can reach the inertial range of 
turbulence. 

^ There is no Ito-Stratonovich ambiguity in this case because the coefficient of Brownian motion does not depend on Q. 
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where L is the maximum extension of the molecule. The force diverges as the molecular extension approaches L; 
thus the extension of the molecule will be finite and smaller than L. The FENE model is appropriate for synthetic 
polymers such as polyacrylamide and polyethyleneoxide; biological macromolecules such as DNA and polyeptides are 
better described by the worm-like chain model (e.g. Larson 2005). In elongational-flow experiments, the flexibility 
parameter b — l? IQ\ usually ranges from 10^ to 10'' (Larson 2005, Shaqfeh 2005). 

The FENE model ([T]) neglects hydrodynamic interactions between the segments of the polymer. Inclusion of 
hydrodynamic interactions makes the coil-stretch transition sharper, but do not change the transition qualitatively 
(e.g. Wiest et al. 1989, Hernandez Cifre and Garcia de la Torre 1998, 1999). Furthermore, the dumbbell model without 
hydrodynamic interactions accurately reproduces the extension-strain curve observed in experiments (Perkins et al. 
1997). The accuracy of the model can be ascribed to the cancellation between the effect of the distribution of drag 
forces along the chain and the effect of the increase in effective drag coefficient with polymer extension (Larson et al. 
1997). Neglecting hydrodynamic interactions has the benefits of providing analytical results. 



III. RELAXATION DYNAMICS 



A steady planar elongational flow is characterised by one direction of stretch, one direction of compression, and one 
neutral direction. The velocity gradient is constant along the directions both of the extensional and compressional 
axes: VjVi = X{SiiSji — 5i2<5j2), A > 0. When a polymer is immersed in that flow the first component of the 
separation vector rapidly becomes much greater than the other components: Qi 3> Q2 and Qi ^ Q3. The extension 
of the molecule, therefore, is approximatively Q ~ Qi, and f{Q) can be replaced by f{Qi).^ According to this 
approximation, the first component of ([1]) gives a one-dimensional stochastic differential equation for the extension 
of the dumbbell: 



dQ = XQ ds 




(2) 



where B{s) is the one-dimensional Brownian motion. For the sake of notational simplicity, we introduce the rescaled 
separation vector q = Q/L. The PDF oi q = \q\, ipil'^t) with q G [0, 1], satisfies the one-dimensional Fokker-Planck 
equation associated with ([2]) (e.g. Stratonovich 1963) 



oq 



Wi - 



qip 



1 d'^i/j 

26 a?"' 



(3) 



where t = s/t, f{q) = f{Lq), and Wi ~ At. The dimensionless number Wi is known as the Weissenberg number 
and measures the level of polymer stretching. For Wi < 1/2 polymers are in the coiled state; for Wi> 1/2 polymers 
are fully extended. The critical value Wic — AcT =1/2 marks the coil-stretch transition in elongational flows (de 
Gennes 1974). Strictly speaking, the approximation leading to ^ and ^ holds true when the polymer is sufficiently 
stretched, i.e. only for Wi ^ Wic. Equation ([3]), indeed, does not yield a good approximation of tp{q;t) at low Wi. 
However, we shall see at the end of this section that it is appropriate to investigate polymer relaxation dynamics in 
terms of © and ^ also below the coil-stretch transition. 

The Fokker-Planck equation ^ is solved with reflecting boundary conditions, that is the probability current 



j{q;t)=Wzqi;{q;t)-l^qi;{q;t)-^^^ 



q,t 



is assumed to vanish at the endpoints of the interval of definition: j{Q;t) 
Under these conditions, the stationary distribution of the extension, ipoi^q) - 
integration (Bird et al. 1987): 



- j{l;t) ^ OWt > (Stratonovich 1963). 
limt^oo -iplq; t), can be derived by simple 



V'o((?)=iVe^^^''(l-g2)l, 0<g<l, 



(4) 



^ The above approximation holds true for a steady uniaxial extensional flow as well: V jVi = X^SnSji — ^iiijijl — SizSj-i/i). Therefore, 
the results henceforth presented are unchanged for this flow. 
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where 

A^ = 2r((6 + 3)/2)/[V^r(V2+l) X iFi(l/2;(6 + 3)/2;6M)] 

and r and iFi denote the Euler Gamma function and the confluent hypergeometric function, respectively (e.g. Erdelyi 
et al. 1953). Large extensions become more and more probable with increasing Wi in accordance with experiments 
(Perkins et al. 1997). Note however that, at low Wi, equation (|4|) does not constitutes a good description of the PDF 
of small extensions due to the approximation behind ([3]). 

The time behaviour of the system depends on the form of the spectrum of the operator ^ with reflecting boundary 
conditions: 



j^{0)= jf.il) = 0, (6) 

where denotes the probability current associated with the eigenfunction Under conditions ([6]), ^ is non- 
positive deflned and symmetric with respect to the scalar product with weighting function l/tpo- Its eigenvalues /i are 
therefore real and non-negative, /i = being associated with the long-time solution ipo (Stratonovich 1963). If ^ has 
a countable spectrum {0, /ii, /j,2, . . . } with /z^ < /ii+i, then T^ei = t j [i\ is the characteristic time needed for i\){fi\ t) to 
attain its long-time form ipoil) when the initial condition ip[q;0) is taken far from equilibrium (Schenzle and Brand 
1979). We note that while r is the relaxation time of the polymer in the absence of flow, T^^i characterises polymer 
relaxation dynamics in the flow. 

To compute Tioi, we need to solve the eigenvalue problem By making the substitutions z = , iIj^{z) — 

(1 — z)iw^{z), we can rewrite ([S]) in the form 

= 0, (7) 



z z — 1 J ^ z{z — 1) 

where w'^ := dw^/dz, 7 = 1/2, 5 — 6/2, a — (l + b — fi/ Wi)/2, (3 — —bWi, and 1/ = b{^— Wi)/2. The above equation 
is a singly confluent Heun equation in the non-symmetrical canonical form (Decarreau et al. 1978a, Decarreau et al. 
1978b, Slavyanov 1995)*. 

Equation ([7|) has two regular singularities at z = and z = 1 and an irregular singularity of Poincare rank 1 
at z = 00. The characteristic exponents at z = are and 1 — 7; at z = 1, they are and 1 — 5. We neglect the 
case of b integer and odd since it does not have physical relevance. Consequently, there are no logarithmic solutions 
and, for |z| < 1, Wf.{z) can be written as W/i(z) = ao<y9o(^) + doz^^^xoiz) with oq, do complex constants and ipo, xo 
analytic functions of z (equivalently of q^) such that (po{0) 7^ and Xo(0) ^ 0. Similarly, for |z — 1| < 1, w^{z) takes 
the form Wf.(z) — ai(pi(z — 1) -f di{l — z^^'^^'^xii^ ~ 1) with ai, di complex constants and ipi, xi analytic functions 
of z such that Lpi[Q) ^ and xi(0) ^ 0. 

It is easily seen that boundary conditions © can be matched only by those solutions of ([7]) that belong to the 
exponent both in z = and z = 1 (that is both do and di have to be zero)'^. The eigenvalue problem defined by ([5]) 
and ^ is therefore mapped into a CTCP on [0, 1] for the confluent Heun equation (O. To compute the eigenvalues /z, 
we adapt the procedure exploited by Svartholm (1939) and Erdelyi (1942, 1944) to ([7]) (see also Slavyanov 1995, 
Slavyanov and Lay 2000). We then expand in series of Jacobi polynomials having the required characteristic 
exponent at z = and z = 1: 



w 



{z) = ^Cn{^i)Un{z) (8) 



n=0 



with 



un{z) = 2F,{-n,u; + n; 7; z) = (-1)- f'f p(-^'^-i)(2z - 1) 

1 (rt + 7) 



* The worm-like chain model, f(q) = 2/3 — l/{(iq) + l/[6(j'(l — q)^], would lead to a second-order linear differential equation with two 
irregular singularities of Poincare rank 1 and 3, respectively. The FENE model in a (^-correlated random flow can be solved in terms of 
a general Heun equation (Martins Afonso and Vincenzi 2005). 

^ We leave aside the case 6 = 2 since it does not have physical relevance. 
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and w = 7+5—1. The function 2F1 denotes the Gauss hypergeometric function; P, 



. p("-7,7-l) 



are the Jacobi polynomials 



of parameters cj — 7 and 7 — 1 and degree n (e.g. Erdelyi et al. 1953). The conditions for the convergence of ([5]) will 
determine the spectrum of ^ . 
By defining the operators 



'Jiw^^ := z{z - 1) 
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z z-lj 
hw^ := -z{z - l)w'^ 



we can rewrite ([7]) in the form 



+ £^2W^ + [pz + cr)Wf_, = 0, 



(9) 



where e = bWi, p = (1 + &) Wi]/2, a — b{Wi — p,)/2. The polynomials it„ satisfy the differential relations (Erdelyi 
et al. 1953) 



and 



with 



^iMn — n {uj + n) Un 
hUn = AnUn+1 + -B,iU„ + C„U„_i 



(10) 
(11) 



n{n + a;)(n + 7) 
(2n + w)(2n + w + 1)' 



Bn — 



n(ri + a;)((5 — 7) 
(2n + uj- l){2n + LU + 1)' 

n(7i + cij)(n + (5 — 1) 



(2n + u;)(2n + uj - 1) ' 
In addition, the following recurrence relation holds (Erdelyi et al. 1953): 

ZUn — AnUji+l + B^Un + CnUn — 1 



(12) 



with coefficients 



{n + 7)(n + u}) 
{2n + uj){2n + uj+l)' 



Br, 



2n{n + LL>) +j{lu- 1) 
(2n + tj-l)(2n + CL) + l)' 

n(n + 6 — 1) 
~ {2n + LU- l){2n + uj) ' 



Inserting expansion ([5]) in ([5]) and exploiting relations (fTU|) . pT|) . (fT^ yield the following three-term recurrence relation 
for c„(^): 



.g_iCo = 
goci + hoco = 

gnCn+l + hnCn + knCn-1 =0 for n > 2, 

where the coefficients 

ffn = sCn+1 + pCn+1, K = n{uj + n) + eBr, + pB,-, + cr, 
depend on fi through p and a. The asymptotic behaviour of hn/gn and kn/gn is given by 

/in 

Sn 



(13) 
(14) 
(15) 



An kn 
and — 

e 9n 



-1 (n — > oo). 



(16) 



According to the Pcrron-Kreuser theorem (e.g. Gautschi 1967, Wimp 1984), implies that the recurrence rela- 
tion (jl5[) has two linearly independent solutions (c+)„ and (c^)n such that 



Cn+i 4n c, 
— ^r- ~ — and 



n+l ^ 

Cn 4n 



(n ^ oo). 



(17) 
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The above sequences have the property: Um„^oo c~/c^ — 0; hence (c~)„ is what is called a minimal solution of (|15p . 
Every other solution non proportional to (c~)„ is asymptotically proportional to (c+)„, and therefore diverges with 
increasing n. Therefore, we can already say that ([5]) can converge only if the sequence (c„)„ is a minimal solution 
of dH). 



The second limit in (|17p implies that lim„^oo \J\cn \ = 0. This condition ensures that, if the coefficients c„ form 
a minimal solution of (jlSp . then expansion ([H]) convergences absolutely to an analytic function in the whole complex 
plane (Szego 1939, p. 252). 

Nonetheless, it is also required that c„ = for n < 0, that is the sequence (c„)„ is subject to initial conditions psp 
and (fT4|) . Equation (fT3|) is trivially satisfied since g-i — 0, and so does (fT4| for /i = 0. (This is in accordance with 
the fact that '00 exists for all Wi and 6.) However, for /i 7^ 0, equation fixes the ratio ci/cq = —(6 + 3)/2. This 
latter requirement is not satisfied by whatever /i, and therefore select the spectrum of ^ . 

To summarise, the spectrum of ^ is the set of those fi such that (c„(/i))„>o is a minimal solution of (fTSj) satisfy- 
ing ci(/i)/co(/^) = —(6 + 3)/2. From Pincherle's theorem (e.g. Gautschi 1967, p. 31), this conclusion is equivalent to 
stating that the eigenvalues ji are the solutions of 

ki/gi 6 + 3 



, / ^2/.g2 

hi/gi 



h2ig2 - 



^3/53 - • ■ • 

The latter represents a transcendental equation for the relaxation spectrum associated with the Fokker-Planck equa- 
tion Q. The reciprocal of the smallest nonzero solution is Troi/r. 

We calculated Trci numerically by means of the modified Miller algorithm (Wimp 1984, pp. 82-85). The behaviour 
of Trei as a function of Wi is reported in fig. [TJ For very small Wi the influence of the external flow is negligible 
and Trei is approximately equal to r. With increasing Wi the ratio T^c\/t starts growing as 1/(1 — 2Wi). This 
behaviour can be deduced by replacing the spring force with a Hookean force, f{q) = 1, and computing the relaxation 
of the moments of q. For a fixed b, the rescaled relaxation time reaches a sharp maximum Tmax in the neighbourhood 
of the coil-stretch transition ( Wic — 1/2); the value of T^ax increases linearly with ^/b (fig. [2]). With increasing 6, Tioi 
attains its maximum value T^ax closer and closer to Wic and the width of the peak decreases (fig. ^ . For large Wi 
the relaxation time is fixed by the time scale of the flow A~^, and therefore Tj-ei/r decreases as Wi^^. 

The computation of T^ei shows that, near the coil-stretch transition, the typical time scale involved in polymer 
relaxation dynamics is strongly different from r. The time associated with the fundamental oscillation mode is not 
representative of polymer temporal dynamics in flow. These results confirm the behaviour encountered in isotropic 
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FIG. 2: (a) Maximum rescaled relaxation time Tma^/r vs. the square root of the flexibility parameter b — /Qq: Tmax/T 
increases as 9i\/b + 62 with 9i ~ 0.69 and 62 ~ —1.18; (6) Behaviour of Wimax as a function of where lyimax is the 
Weissenberg number at which the relaxation time attains its maximum value Tmax; (c) Relaxation time vs. Wi for different b, 
rescaled to make the maximum values coincide. 



(5-correlated Gaussian flows (Celani et al. 2005a, Martins Afonso and Vincenzi 2005), but here the enhancement of Tioi 
is stronger since the elongational flow is more efficient in deforming polymers. 

In Brownian Dynamics simulations and experiments the relaxation time of 'ip[q]t) can be measured from the decay 
of the moments of the extension to their steady-state value. To observe Trci, one has to make sure that the initial 
condition il){q] 0) is not orthogonal to the first eigenfunction of .Sf, that is ■i/'i (9) in ©, with respect to the scalar product 
with weighting function ilJo{q). Figure [3] shows the shape of i^iiq) for different Wi. An initial PDF concentrated either 
at a coiled or a stretched extension q (formally ipil'^O) = ^{q ~ q)) has projection onto ^i{q) equal to "01 (g)/i/'o (9) 
by construction (e.g. Stratonovich 1963). The initial condition with all the polymers having the same extension q is 
therefore not orthogonal to ipi{q) if ipi{q) ^ 0; all the moments will thus tend to their stationary value with a typical 
time scale Trei. 

We now show that it is accurate to compute Trd by means of ([3]) also for Wi < Wic- At low Wi, the entropic force 
can be modelled as a Hookean force, i.e. f{Q) = 1 in (jT]). The vector equation ([T|) therefore reduces to a set of three 
decoupled stochastic differential equations for Qi, Q2, Q3. Hence we have ^{q;t) = ^p^^'' {qi;t)^p^'^\q2;t)^p^^\q3;t), 
where \I'(q;t) is the PDF of the vector q and ^'''^\qi;t) is the PDF of the i-th component of q. The relaxation 
time of ^{q;t) is thus the longest of the relaxation times of ijj^^\qi; t) , ii)^'^\q2]t), and ip^^^q^i^t) given that all 
the ^p^'^'' {qi;t) have a non-trivial long-time limit. It is easy to check that at low Wi the relaxation time of ip^^\qi;t) 

is = t/(1 - 2 Wi), while T^^^ = r/(l + 2 Wi) and T^^^ = t (with obvious notation). The function ■^{q; t) therefore 
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FIG. 3: Stationary PDF of the extension and eigenfunction of the Fokker-Planck operator associated with the smallest 
nonzero eigenvalue [i\ {h = 600.25). The eigenfunctions have been computed numerically by means of the iteration-variation 
method (Morse and Feshbach 1953). 



has the same relaxation time as ■i/'^^-' Now we note that ipiQ'jt) a-nd have the same long-time behaviour 

since ip{q;t) = J^^ 'i'{q{q,d, <j))',t) sin-d dddcf). We thus come to the following conclusion: the relaxation times 

ofip{q',t) and ip'^^Hlu coincide at low Wi and t/j^-^'^ {qi; t) satisfies ([3]). This fact explains why ([3]) yields an accurate 
computation of Trci also below the coil-stretch transition. 



IV. SUMMARY AND CONCLUSIONS 



We have investigated the relaxation dynamics of an isolated polymer in an external flow. Previous studies focused 
on the relaxation of a polymer in a solvent driven only by Brownian motion. Those studies determined the time t 
associated with the slowest oscillation mode of the molecule in the absence of external flow. We have considered an 
elongational flow and derived a transcendental equation for the relaxation spectrum associated with the time evolution 
of the PDF of polymer extension. The problem of computing this spectrum has been recast as a central two-point 
connection problem for a confluent Heun equation. The Heun equation results from the separation of the time variable 
in the Fokker-Planck equation for the PDF of the extension. 

Our analysis shows that an external elongational flow strongly influences polymer relaxation dynamics. The longest 
relaxation time associated with the evolution of the PDF of the extension provides an estimation of the time scale of 
polymer deformation. Near the coil-stretch transition, this time is significantly greater than t (already one order of 
magnitude greater for short molecules). The physical reason is the large number of polymer configurations close to the 
coil-stretch transition. At intermediate Weissenberg numbers polymer conformation results from a critical competition 
between the entropic force and the velocity gradient. Therefore, around the transition, coiled and stretched polymer 
coexist in the flow, and this makes the relaxation to the equilibrium PDF particularly long. It is worth noticing 
that the ampliflcation of the relaxation time is stronger when hydrodynamic interactions are taken into account; this 
behaviour is intimately related to the finite-time conformation hysteresis observed in elongational fiows (Celani et al. 
2006). 
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